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Abstract. For future configurations, we study the relation between the 
. abatement of the noise sources and the Signal to Noise Ratio (SNR) for coalescing 

I _ ' binaries. Our aim is not the proposition of a new design, but an indication of 

>«-' ' where in the bandwidth or for which noise source, a noise reduction would be 

>D ' most eflicient. We take VIRGO as the reference for our considerations, solely 

(_P ^ applicable to the inspiralling phase of a coalescing binary. Thus, only neutron 

^^ , stars and small black holes of few solar masses are encompassed by our analysis. 

C 3 ■ The contributions to the SNR given by final merge and quasi-normal ringing are 

neglected. It is identified that i) the reduction in the mirror thermal noise band 
provides the highest gain for the SNR, when the VIRGO bandwidth is divided 
according to the dominant noises; ii) it exists a specific frequency at which lies 
^^ , the potential largest increment in the SNR, and that the enlargement of the 

{3J[). bandwidth, where the noise is reduced, produces a shift of such optimal frequency 

to higher values; iii) the abatement of the pendulum thermal noise provides the 

largest, but modest, gain, when noise sources are considered separately. Our 

K^ recent astrophysical analysis on event rates for neutron stars leads to a detection 

KN . rate of one every 148 or 125 years for VIRGO and LIGO, respectively, while 

a recently proposed and improved, but still conservative, VIRGO configuration 
would provide an increase to 1.5 events per year. Instead, a bi-monthly event 
rate, similar to advanced LIGO, requires a 16 times gain. We analyse the 3D 
(pendulum, mirror, shot noises) parameter space showing how such gain could be 
achieved. 



1. Introduction 

Improved detectors are currently under consideration in Europe. The motivation 
of this work is the analysis of the cost/benefit ratio in support to the forthcoming 
R&D effort, a major concern being how to gear efforts such that the improvements 
will sensibly increase the probability of detection and give birth to an observational 
relativistic astrophysics. To such end, this work takes a very pragmatic approach. The 
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design of future interferometers is well beyond the scope of this paper and instead we 
deal with the detector from an user view point. Once the most efficient utilisation 
strategy shall be identified, detailed design and R&D should follow accordingly. In this 
respect, the approach adopted herein proposes the end-user, i.e. the astrophysicist, 
to guide the requirements on the design of the interferometers. A not very different 
approach was taken by Babusci and Giovannini (2001) who investigated stochastic 
and isotropic gravitational wave backgrounds produced in pre-big-bang models. 
Scientific interest, knowledge of the model, and number of events hint to coalescing 
binaries as primary candidates. Thus, improvements in signal to noise ratio (SNR) by 
noise abatement, for inspiralling binaries of small chirp mass, are discussed. For this 
exercise, we make reference to the VIRGO interferometer. 

The structure of the paper is the following. In the next section, we show the event 
rates out of a recent analysis (de Freitas Pacheco et al., 2004). In the third section, 
we briefly summarise the main expressions that determine the SNR for inspiralling 
binaries and present the main features of the VIRGO sensitivity curve. In the fourth 
section, we perform the analysis leading to our main findings on the relation between 
noise abatement and the increment in the SNR. In the fifth section, we estimate the 
detection rate of a recently proposed VIRGO configuration and address the necessary 
improvements to acquire a detection rate similar to advanced LIGO. 

2. Event rates 

The estimates of the coalescence rate of NS-NS binaries are performed in two main 
steps. Firstly, the merging rate in our Galaxy is evaluated and then, by assuming 
that this value is typical, extrapolations to the local universe are performed under 
the assumption of some adequate scaling. The galactic rate has been estimated by 
different authors and ranges from 10~^ /yr~^ up to few times 10~^ /yr~^, but values 
as high as (2 — 3) x IQ^^ /yr~^ have already been reported (Tutukov & Yungelson 
1993; Lipunov 1997; Kalogera et al. 2004). Here we adopt the approach by de Freitas 
Pacheco (1997) and de Freitas Pacheco et al. (2004). 

Let us suppose a massive binary (masses of the components higher than 9 Mq) 
formed at the instant t' . Let r* be the mean evolutionary timescale for the system 
to evolve into two neutron stars, typically of the order of 10* yr. Define P{t) as 
the probability per unit of time for a newly formed NS-NS binary to coalesce in a 
timescale r and define i?*(t) as the star formation rate (in Moyr"^). Then, the 
galactic coalescence rate at instant t is 

/•(t-Tj-To) 

Vc{t)^Mns\ P{T)R,{t - T, - Tn)dT (1) 



where fb is the fraction of massive binaries formed among all stars, (3ns is the fraction 
of formed binaries which remain bounded after the second supernova event and 
A is the fraction of formed stars in the mass interval 9-40 Mq. We have set tq 
as the minimum timescale for a NS-NS binary to coalesce. Numerical simulations 
discussed in de Freitas Pacheco et al. (2004) allow the determination of tq, the 
product fbPns and the probability distribution P{t). One obtains Tq = 2 x 10^ yr, 
P{t) oc 1/t and fbPns — 0.00326. This latter value corresponds to a natal kick velocity 
dispersion of about 80 km/s. If we adopt the galactic star formation history derived 
by Rocha-Pinto et al. (2000), it results for the present galactic NS-NS coalescence rate 
fs = (1-7 ± 1.0) X 10~^ yr~^, where the estimated uncertainty is essentially due to 
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uncertainties in the natal kick velocity dispersion and in the actual ratio between the 
number of single to binary pulsars. A similar calculation for a typical elliptical galaxy 
with an absolute magnitude Mb = —20.7 gives a coalescence rate i^e = 8.6 x 10~^ yr~^ 
(de Freitas Pacheco et al. 2004). The weighted local coalescence rate was obtained 
assuming a fraction of 65% of spirals and 35% of ellipticals and lenticulars (SO's) and an 
elliptical-to-spiral luminosity ratio of 1.26. The resulting rate is i^c = 3.4 x 10~^ yr~^. 
Extrapolation from the local coalescence rate to the expected rate within a 
spherical volume of radius D probed by the detector is made through the scale factor 
Kb{D). This factor is defined as the ratio between the total blue luminosity within 
the considered volume and the Milk Way luminosity. We use here the scale factor 
computed by de Freitas Pacheco et al. (2004), who considered the distribution of 
galaxies derived from the LEDA databasis and included the contribution of the Great 
Attractor, centered at the Norma cluster. The expected rates as a function of distance 
are given in tab. 1. 



D (Mpc) 


Event rate yr ^ 


10 


3.44 X 10-^ 


15 


2.00 X 10-^ 


20 


3.00 X 10-^ 


30 


6.44 X 10-^ 


40 


1.21 X 10"^ 


50 


2.20 X 10-^ 


70 


7.08 X 10-^ 


100 


1.63 


150 


2.30 


200 


5.44 


300 


18.4 



Tab. 1 Expected event rates as a function of distance. 



3. Inspiralling binaries and VIRGO features 

Matched filtering is the baseline detection tool, although computationally demanding. 
It enhances the signal by the square root of the number of observed cycles in the binary. 
For our analysis, the simplest post-Newtonian model is assumed. A more accurate 
model may marginally change the absolute values of the SNR but not their relative 
values we are interested in. A simplified expression, derived by the Peters-Mathews 
model (Peters, 1964; and Mathews, 1963) is given by (Spallicci et al., 1997): 




1 



-df 



where /(/) is the one-sided spectral density of the noise, 
T 1 TW^/e 



A 



M - (/iM2/3)3/5 



(MiMa)^/^ 
(Ml + M2)i/5 



(2) 



(3) 



M,iJ,,M, Ml, M2 being the chirp, reduced, total and single masses of the binary. The 
coefficient cr = 2/5 takes into account the average relative orientation between binary 
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and detector (the discrepancy with more accurate expressions is around 5%). The 
VIRGO sensitivity curve, fig. 1, is given by§: 

loidif) = 3.8 X lO^^V^ + 9.5 X 10""*^ + 1.6 x 10"''^"^ + 8.3 x lO^^V"^ (4) 

and it is divided into the foUowing zones dominated by a specific noise: a. 1 - 2.3 Hz 
seismic, b. 2.3 - 52 Hz thermal pendulum, c. 52 - 148 Hz thermal mirror, d. 148 - 
10.000 Hz shot. We define (Spallicci, 2003) the gain simply as (the upper limit of the 
integrals is known to be irrelevant if sufficiently large to encompass the total inspiral): 



^ _ {SNR),,, 



{SNR)oid 



I f'° ^^- df N^ '^' 

1 Hz J^newJ 

10 KH. ^j 



(5) 



V A H. loMp/^ I 

,(/) generally indicates an improved spectral noise density. 

4. Analysis for an improved SNR 

We emphasize that eq. (jSj) is a ratio and that most of the consequences of simplifying 
assumptions (on the antenna and the source) will be not influential. Our considerations 
are applicable solely to the inspiralling phase of a coalescing binary. This is equivalent 
to say that only neutron stars and small black holes of few solar masses are 
encompassed by this analysis. The contributions to the SNR given by final merge 
and quasi-normal ringing are neglected. 
We introduce the noise abatement factor fc^ for which: 



1 



loid kj 



[u{f + F)-u if -F)] + [l-uif + F)+u if - F)] (6) 



where 



F^^-fo uif) = 0(/ < 0), 0.5(/ = 0), 1(/ > 0) 

We note that the Gate function is built on the Heaviside Step u{f) function. The 
operator kf acts on a sub-band A/ centered on /q. Therefore, the function Inew takes 
the values loid or loid/k'i, outside and within the chosen window, respectively. 
The first investigation aims to identify the SNR improvement for abatement of the 
noise in one of the four a-d sub-bands defined in the previous section. We improve 
each of the four sub-bands in turns; fig. 2 shows the behavior of the SNR gain for 
different values of noise abatement. We observe that the abatement in the thermal 
mirror noise band (c) provides the largest SNR gain. 

The second investigation aims to identify the optimal frequency around which the 
noise abatement should be performed. We firstly look for the single frequency at 
which the noise reduction is more beneficial, identified as 52 Hz [fig. 3). Further, we 
analyse where the noise reduction in a given bandwidth of 10 or 100 Hz would be more 
beneficial ifig. 3). The increment of the bandwidth shifts the optimal frequency to 
larger values: for a 100 Hz bandwidth, the optimal frequency is 75 Hz. Furthermore, 
equivalent situations achieve the same gain in the SNR offering alternative strategies. 
Fig. 3 reports a similar gain G = 2.5 for fc/ = 3 at 75 Hz in a bandwidth of 100 Hz 

§ http://www.virgo.infn.it 
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or for kf — 10 at 30 Hz in a bandwidth of 10 Hz. The maximum gain is G = 8 for a 

bandwidth of 100 Hz centred at 75 Hz. 

In the 52-148 Hz band, the major hmiting factor is coming from the thermal noise of 

the mirrors but a large noise reduction would imply lowering the contributions from 

the pendulum and the shot noises too, since a large reduction in thermal noise is 

useless without also reducing the other noise sources. 

If instead, we consider the noise abatement per noise source, the new spectral noise 

density is given by: 



E^: 



hr^, , ^^^ 



having attributed a noise reduction factor fc„ to a single source, out of r + 1 noises. 
Indeed, the sole reduction of the thermal mirror noise leads to a gain G — 1.35, fig. 4- 
When we consider the four major noise sources and apply to each of them, in turns, 
a reduction factor A:„, fig. 4-, the gain G is never larger than 1.55 (thermal pendulum 
noise). Obviously, asymptotic gain values are achieved, whenever only one single noise 
source is dealt with, while the others are left unchanged. 

5. Detection rates 

The actual specifications on the VIRGO sensitivity show that a SNR of 7 is achieved 
at 13 Mpc for neutron stars binaries combined with a false alarm rate (FAR) of 1 per 
year. This distance corresponds to an event rate of 1 every 148 years (tab. 1), while 
the event rate for LIGO is 1 each 125 years (Pacheco et al., 2004), same SNR and 
FAR being at 14 Mpc. A gain G of about 8 in the SNR, fig. 3 would allow VIRGO to 
explore the universe up to about 100 Mpc with the same SNR and FAR. To the latter 
values, a detection rate of 1.5 event per year is associated. 

Recently, an improved configuration for VIRGO has been proposed (Punturo, 2004), 
its sensitivity curve being shown in fig. 5. The new sensitivity is obtained by reducing 
the pendulum noise by a factor 28, the thermal noise by 7 and the shot noise by 
4, throughout the spectrum. The associated gain G is 7.5. The Punturo (2004) 
configuration is rather conservative and is currently analysed for further improvement. 
Future configurations will likely consider a cryogenic design||. From tab. 1, we evince 

II The suspended test masses thermal noise power spectrum depends upon the temperature, the strain 
energy and the loss angle of the mirror (Harry et al., 2002; Bondu et al., 1998). As demonstrated 
by Braginsky et al. (1998, 2003), the thermo-elastic effect is present both in the mirror and in 
the coating (thermo-dynamical contribution) . The direct temperature dependence of the fluctuation 
power spectrum suggests that, if we want gain an order of magnitude in terms of the space-time 
strain h, we must decrease the temperature of about two orders of magnitude, that is go to cryogenic 
temperatures. Planning entails use of crystalline Silicon to realize the suspension for its lower loss 
angle than fused silica; its thermal expansion coefficient goes essentially to zero before 10 K. Also, 
the thermal conductivity of the Silicon increases at low temperature pushing up in frequency the 
residual thermo-elastic peak. The latter property of Si is also fundamental for its use in an advanced 
interforometric detector, in which the power circulating is so high that the thermal lensing effect 
in the mirrors could change the optical properties of the Fabry-Perot cavities, affecting the whole 
functioning. The Silicon thermo-mechanical properties suggest it as a good candidate also for the 
mirrors. All the thermo-dynamical phenomena related to the substrate would disappear using Silicon 
mirrors at cryogenic temperatures. Different is the behavior of the coatings. Recent preliminary 
measurements showed that the loss angle of Ta2 05 coatings doesn't decrease with the temperature. 
In any case, at least a reduction of a factor 10 due to the direct temperature dependence is expected 
in terms of noise amplitude spectral density. Obviously, in order to use Silicon mirrors, it is necessary 
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that it is mandatory that a detector reaches at least 80 Mpc to get an yearly rate 
of events. A significant further relevant improvement imposes VIRGO to go beyond 
210 Mpc for one event each two months, as advanced LIGO. Fig. 6 shows which 
combinations of concurring abatements on the three most dominant noises (thermal 
pendulum, thermal mirror, shot) produces a gain G = 16, necessary to achieve 210 
Mpc with SNR = 7 with an yearly FAR. The dots represent integer values of the noise 
reduction factors for thermal pendulum and mirror, joined by isolines of the shot noise 
reduction factor. 

6. Conclusions 

We have analysed the efficiency of potential improvements on the actual VIRGO 
configuration versus the SNR for coalescing binaries. Recently proposed configurations 
would explore a sphere of less than 100 Mpc of radius (SNR — 7, yearly FAR) and 
provide 1.5 detections per year, while the actual configuration would detect 1 event 
each 148 years, exploring up to 13 Mpc (SNR — 7, yearly FAR). A bi-monthly rate 
requires an improvement of a factor 16 in gain. The noise reduction parameter space 
necessary to achieve such gain has been explored. 
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Figure 1. VIRGO strain sensitivity h{f)/\^Hz. 
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Figure 2. Gain G versus noise reduction factor kf for the four sub- bands a-d 
(a. 1 - 2.3 Hz seismic, b. 2.3 - 52 Hz thermal pendulum, c. 52 - 148 Hz thermal 
mirror, d. 148 - 10.000 Hz shot). 
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Figure 3. Gain G in the frequency range 1-250 Hz for three levels of abatement 
of noise {kf = 3,5, 10), relative to an improvement in a bandwidth of 1, 10, 100 
Hz. 
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Figure 4. Gain G versus abatement of noise fc„ for tlie four major noises 
(Newtonian, thermal pendulum, thermal mirror, shot). 
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Figure 5. VIRGO strain sensitivities h{f)/\/Hz: the actual configuration and 
an example of an improved configuration (Punturo, 2004). 
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Figure 6. The dots correspond to the integer values of the kn reduction factors 
for thermal pendulum and mirror noises, parametrised by the shot noise. Any 
threefold combination leads to an equal gain G = 16 for VIRGO and therefore a 
bi-monthly detection rate for NS-NS coalescences, as advanced LIGO. 



